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Coulomb corrections to perturbation theory for producing electron- 
positron pairs in ultrarelativistic heavy ion collisions are considered in a part- 
analytical, part-numerical approach. Production probabilities are reduced 
from perturbation theory with increasing charge of the colliding heavy ions, 
as has been previously argued in the literature. It is shown here that the re- 
duction from perturbation theory comes from the appropriate physical spatial 
cutoff of the electromagnetic potentials arising from the colliding ultrarela- 
tivistic heavy ions. 
PACS: 25.75.-q, 34.90.+q 

I. INTRODUCTION 

The problem of calculating heavy ion induced continuum e + e~ pair production to all 
orders in Za has received some renewed interest in the past several years. Realization that 
in an appropriate gauge [1], the electromagnetic field of a relativistic heavy ion is to a very 
good approximation a delta function in the direction of motion of the heavy ion times the two 
dimensional solution of Maxwell's equations in the transverse direction [2], led to an exact 
solution of the appropriate Dirac equation for excitation of bound-electron positron pairs 
[3]. Given this solution, it was perhaps not surprising that the solution of the Dirac equation 
was obtained independently and practically simultaneously by two different collaborations 
[4-6] for the analagous case of continuum e + e~ pair production induced by the corresponding 
countermoving delta function potentials produced by ultrarelativistic heavy ions in a collider 
such as RHIC. An extended discussion and reanalyis of this solution, with comments on early 
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parallel work in the literature, shortly followed [7]. One apparent physical consequence of 
this solution was that the rates for pair production in the exact solution agreed with the 
corresponding perturbation theory result [5-7]. 

Several authors subsequently argued [8-10] that a correct regularization of the exact 
Dirac equation amplitude should lead to deviations from perturbation theory, the so called 
Coulomb corrections. Although, as has been pointed out [11], the derived exact semiclas- 
sical Dirac amplitude is not simply the exact amplitude for the excitation of a particular 
(correlated) electron-positron pair, there are observables, such as the total pair production 
cross section, that can be constructed from this derived amplitude. The exact amplitude for 
a correlated electron-positron pair will not be treated here. It is the Coulomb corrections to 
the observables that can be constructed from this exact Dirac equation amplitude that are 
the topic of this paper. 

In what follows it will be shown from a somewhat different approach from what has been 
done before that Coulomb corrections must exist, that they arise from the physical cutoff of 
the tranverse Coulomb potential, and the accuracy of their evaluation has been up to now 
limited by an effective two-peak approximation to the exact retarded Dirac amplitude. 



One begins the semiclassical Dirac solution by representing the electromagnetic effect of 
one heavy ion on the other as the Lienhard-Wiechart potential produced by a point charge 
on a straight-line trajectory 



b is the impact parameter, perpendicular to the z-axis along which the ions travel, p, z, 
and t are the coordinates of the potential relative to a fixed target (or ion), a z is the Dirac 
matrix, and Z, v and 7 are the charge, velocity and relativistic 7 factor of the moving ion. 
If one makes a gauge transformation on the wave function [1] 



II. THE DIRAC EQUATION SOLUTION 




aZ(l — va z ) 



(1) 
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where 



aZ 



In [7 (z - vt) + ^b 2 + 1 2 (z-vt) 2 } 



(3) 



v 



the interaction potential V(p,z,t) is gauge transformed to 



V(p,;M) 



aZ{l — va z ) aZ(l — (l/v)a z ) 



(4) 



The second term is pure gauge and serves to reduce the range of the potential in (z — vt) 
to more closely map the (z — vt) range of the B and E fields, which have the denomenator 



This is the potential that allowed the closed form solution of the Dirac equation for the 
bound-electron positron problem. The full solution of the problem is in perturbation theory- 
form, but with an eikonalized interaction in the transverse direction 



in place of the pertubation interaction Eq. (5) producing the higher order effect in Za. 
Recall that this exact semiclassical solution produced a reduction of a little less than 10% 
in the predicted cross section for Au + Au at RHIC [3]; one can identify this reduction as 
a Coulomb correction to bound-electron positron pair production. 

In the bound-electron positron problem one conveniently takes the electromagnetic field 
of one moving heavy ion seen in the rest frame of the heavy ion that receives the created 
electron. For production of continuum pairs in an ultrarelativistic heavy ion reaction one 
may work in in the center of mass frame and the electromagnetic interaction goes to the 
limit of two countermoving 8 function potentials 




V(p, z, t) = -5{z - t) (1 - a z )aZ P In (b - pf . 



(5) 



V(p,z,t) = -iS(z - t)(l - a z )(exp[-iaZ P \n(b - pf] - 1). 



(6) 
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V(p, z, t) = 5{z -t)(l- a z )A-(p) + 8(z + t)(l + a z )A + (p) 



(7) 



where 



(p) = — Za In 



(p±b/2) 2 

(6/2) 2 



(8) 



The semi-classical Dirac equation with this potential has been solved in closed form [4-7] . 
Baltz and McLerran [5] noted the apparent agreement of the obtained amplitude with that 
of perturbation theory even for large Z. Segev and Wells [6] also noted the agreement with 
perturbation theory and noted the scaling with Z\Z\ seen in CERN SPS data [12]. These 
data were obtained from reactions of 160 GeV/c Pb ions on C, Al, Pa, and Au targets as well 
as 200 Gev/c S ions on the same C, Al, Pa, and Au targets. The group presenting the CERN 
data, Vane et al., stated their findings in summary: "Cross sections scale as the product 
of the squares of the projectile and target nuclear charges." On the other hand, it is well 
known that photoproduction of e + e~ pairs on a heavy target shows a negative (Coulomb) 
correction proportional to Z 2 that is well described by the Bethe-Maximon theory [13]. 



As noted in the Introduction, several authors have argued that a correct regularization of 
the exact Dirac equation amplitude must lead to Coulomb corrections. The first analysis was 
done in a Weizsacker- Williams approximation [8]. Subsequently, Lee and Milstein argued 
[9,10] the existence of Coulomb corrections by an approximate analysis of the closed form 
solution of the Dirac equation. We will take as our starting point a somewhat extended 
consideration of the results of Lee and Milstein. 

To begin let us write the previously derived semiclassical amplitude for electron-positron 
pair production [4-7] in the notation of Lee and Milstein [9] 



p and q are the four-momenta of the produced electron and positron respectively, k is an 
intermediate transverse photon momentum to be integrated over, 
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III. COULOMB CORRECTIONS 





ka<\ \ -i \ a(k- p ± ) +7 m 

M[k) = u(p) - j-— 7_w(-g) 

—p+q~ — (k — p^) 2 — m 2 + «e 

+w(p) 7j vj ^— — 7+«(-5) (10) 

—p-q+ — (k — qjj 2 — m 2 + 

and the effect of the potential Eq. (7-8) is contained in integrals, F B and Fa, over the 
transverse spatial coordinates taking the form 

F(k) = J d 2 p exp[-ik • p]{exp[~i2Zalnp] - 1} 

= 2tt / pdpJ (kp){exp[-i2Za\np] - 1}. (11) 

JO 

F(k) has to be regularized or cut off at large p. How it is regularized is the key to under- 
standing Coulomb corrections. If one merely regularizes the integral itself at large p one 
obtains [5-7] apart from a trivial phase 



AiraZ 

F M = p^z (12) 



All the higher order Za effects in M(p, q) are contained only in the phase of the denom- 
inator of Eq. (12). As we will see, it directly follows that calculable observables are equal 
to perturbative results. 



A. Observables 

Before considering the Lee and Milstein analysis, we will discuss the observables that 
can be calculated [14-17] from the solution of a Dirac equation such as Eq. (9-12). We have 
pointed out that the derived semiclassical Dirac amplitude M(p, q) is not simply the exact 
amplitude for the excitation of an electron-positron pair [11]. The point is that exact solution 
of the semi-classical Dirac equation may be used to compute the inclusive average number 
of pairs — not an exclusive amplitude for a particular pair. Calculating the exact exclusive 
amplitude to all orders in Za is not easily tractable due to need for Feynman propagators 
[11]. The possibility of solutions of the semi-classical Dirac equation is connected to the 
retarded propagators involved. In this paper we do not consider the exclusive (Feynman 
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propagator) amplitude at all. We concentrate on observables that can be construted from 
the above amplitude and investigate the Coulomb corrections contained in them. 

The occupation number or inclusive number of electrons created in state p (at impact 
parameter b) is 

N ^ = I H$7, m "' q)l2 (13) 

Likewise the inclusive number of positrons created in state q is 

N ^-I0k mp - q)f (14) 

These inclusive expressions say nothing about correlations between electrons in state p and 
and positrons in state q. 

The mean number of electron-positron pairs is of course equal to either the mean number 
of positrons or the mean number of electrons and may be obtained by integrating over either 
of the previous expressions. 



It is possible to calculate well-defined observables from the occupation numbers by inte- 
grating over the impact parameter b 

da(p) = J d 2 bN(p) = j d 2 b0^\M{p,q)\\ (17) 
da(q) = j d%N(q) = j d%^^\M(p,q)\ 2 . (18) 



and 



o* = /^ = /A$$g\M M f>. (19, 

da(p) is the cross section for an electron of momentum (p) where the state of the positron 
is unspecified. Likewise da(q) is the cross section for a positron of momentum (g) with the 
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state of the electron unspecified. Note that ot corresponds to a peculiar type of inclusive 
cross section which we should call the "number weighted total cross section" , 

oo 

d 2 bN= / rf 2 6^nP n (6), (20) 

J n=l 

in contrast to the usual definition of an inclusive total cross section 07 for pair production, 

„ 00 

a I = j d 2 bY,Pn(b). (21) 



n=l 



Now we can write for the factor common to all the cross sections 

/^|M(p, 9 )r = /^/|l/gexp[ i (k-k0.b] 
x M(k)M(k'YF B (k)F B (k'y 

X F A (q ± + p ± - k)F A (q ± + p ± - k')*. (22) 

Integrating exp[i(k — k') • b] over the impact parameter b in the usual way gives 
(2vr) 2 5(k - k') and so 

/ d 2 b\M(p,q)\ 2 = J ^\M(k)\ 2 \F A ( q± + P± -k)\ 2 \F B (k)\ 2 . (23) 

One now obtains expressions for da(p), da(q), and <jt that appear identical to the result 
of perturbation theory (scaling as Z\Z B ) when our previous expression for F(k) Eq. (12) 
is employed. 



= / Wkj + Px - k)| V S (k)! 2 (25, 

Obviously F B and F A still have to be regularized or cut off at small |k| and |qj_ + pj_ — k|. 
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B. The regularization of Lee and Milstein 



The strategy of the first paper of Lee and Milstein [9] was to evaluate Coulomb corrections 
by Taylor expanding M around k = 0, i.e. M(k) ~ k • L. The derivative L is evaluated 
at k = 0, and also in the evaluation of e.g. Eq. (26) k is ignored in F^(q± + pj_ — k). All 
the k dependence of the integral is then contained in d 2 k k 2 \F B (k)\ 2 . Lee and Milstein then 
invite us to consider the integral representing the difference between the exact solution and 
the perturbative solution 




where 

F(k) = J d 2 p expf-ik • p]{expH X (p)] - 1}, (28) 
with the transverse form of the potential not yet specified 

/oo 
dzV(z,p), (29) 
-oo 

and 

F°(k) = -ij d 2 p expHk • p] X (p) (30) 

is the perturbative expression limit of F(k). 

Lee and Milstein keep the 2Za\n(p) form for x{p) but switch the order of integration 
between p and k. They integrate k to some finite upper limit Q an then claim to set Q 
to infinity in the resulting expression. Actually Q simply falls out of the problem by a 
rescaling of p to p/Q. Next, after integrating over the rescaled p, the expression they obtain 
is a universal function of Za 

G = -8n(Za) 2 [Re^((l + iZa) + lEu ierl (31) 

where ^((l + iZa) is the digamma function and ^Euier is Euler's constant. This expression 
may be alternatively expressed as 
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G = -8n(Za) 2 f(Za), (32) 

where f(Za) is the same function that was derived by Bethe and Maximon for Coulomb 
corrections to e + e~ photoproduction on heavy nuclei and takes the form 

OO 1 

f(Za) = (Za) 2 £ (33) 

The derivation and result may seem a little mystifying. Lee and Milstein state, "Thus, 
we come to a remarkable conclusion: although the main contribution to the integral in Eq. 
(4) comes from the region of small k, where |F(k)| differs from (|F°(k|) = AnZa/k 2 and 
depends on the regularization parameters (the radius of screening), nevertheless the integral 
G itself is a universal function of Za." As we will see later the only part of this quoted 
statement that is completely true is that G is a universal function of Za. 

G is then used by Lee and Milstein to calculate the Coulomb correction arising from ion 
B by taking ion A is to lowest order in Za. Generalizing this approach, the corresponding 
Coulomb correction arising from ion A is also evaluated [10]. The sum of these two contri- 
butions then agree with the Coulomb corrections as evaluated by Ivanov, Schiller, and Serbo 
[8] using the Weizsacker- Williams method. 



C. A physical regularization 

Let us try to understand Lee and Milstein's result by putting in a physical cutoff to the 
transverse potential x(p) (which has been up to now set to 2Za In p). Instead of regularizing 
the integral itself and letting the cutoff radius go to infinity as was originally done [4-7], 
we will apply an appropriate physical cutoff to the interaction potential. In the Weizsacker- 
Williams or equivalent photon treatment of electromagnetic interactions the potential is cut 
off at impact parameter b ~ 7/0;, where 7 is the relativistic boost of the ion producing the 
photon and uo is the energy of the photon. As Lee and Milstein subsequently recall (but do 
not utilize) if 

X (p) = r dzV^z 2 + p 2 ) (34) 

J— CO 
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and V(r) is cut off in a physically motivated way, such as an equivalent photon cutoff, then 

V(r) = -^M-r^s/l] (35) 

where 

p+ + <?+ p- + q- 

"a = 2 ; Wb = 2 ^ ) 

with uja the energy of the photon from ion A moving in the positive z direction and oj b the 
energy of the photon from ion B moving in the negative z direction. For simplicity we will 
suppress the subscripts on u, remembering however for future possible use that uja,b are 
well defined in terms of p± and q±. The integral Eq. (34) can be carried out to obtain 

X (p) = -2ZaK (pu/>y), (37) 

and 

F A)B (k) = 2tt | dppJ (kp){exp[2tZ AB aK (puj/ 1 )} - 1}. (38) 

The modified Bessel function K (pu/^) = — \n(pu/2^) for small p and cuts off exponentially 
at p ~ 7/c<j. This is the physical cutoff to the transverse potential. 
One may define £ = kp and rewrite Eq.(38) 

F a ,bM = %\ dttU0{evp[2iZ A , B aK Q {&hk)] - 1}. (39) 

It is now clear that Fa,b is a function of 1/k 2 times some function of {^fk/u). The pertur- 
bative limit F° B (k) is analytically soluable and takes the form 

= 47iZ A<B a = 47iZ AtB a 

A ' By ' fc 2 + w 2 / 7 2 A; 2 (1+cj 2 /P 7 2 ) 1 J 

Fig. (1) displays the results of numerical calculation of the scaled magnitude of F(k) as 
a function of Ary/cj for Z = 1 (essentially the perturbative form Eq. (40)) and for Z = 82. 
Note that the upper cutoff of p at 7/0; has the effect of regularizing F(k) at small k. F(k) 
goes to the constant 47r 7 2 /o; 2 as k goes to zero in the Z — 1 perturbative case; it goes to a 
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reduced constant value as k goes to zero for Z = 82. The form of the original solution Eq. 
(U) 



F(k) 



AnaZ 



(41) 



fc2-2iaZ 

is simply wrong because it is unphysical. Since it lacks a proper physical cutoff in p, it not 
only blows up at k = 0, but it also fails to exhibit the correct reduction in magnitude that 
occurs when kj/u is not too large. 

Coulomb Corrections 

Scaled Magnitude of Transverse Integral F(k) 

1 



0.5 



0.01 




100 



FIG. 1. The decrease in the magnitude of the transverse integral F with Z. The two sets of 
curves have been normalized to display that the finite Coulomb correction only rescales down the 
|.F(k)| ~ 1/uj 2 behavior at kuo/j = and that the negative Coulomb corrections do not vanish 
until well above the onset of |-F(k)| ~ 1/k 2 dominant behavior. 

Fig. (2) is an alternate display of results of the numerical calculations showing the 
fractional decrease in the ratio |F(k)|/|F°(k)| for various values of Z as a function of kj/u. 
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It is clear from the two Figures that for increasing Z Coulomb corrections reduce F(k) from 
the perturbative result for k^/u « 100. Only for k > ~ 100 cj/7 does the magnitude of 
F(k) go over into the original form of Eq. (41). 

Coulomb Corrections 

Fractional Decrease of IFI from with Z 

0.15 | 




ky/co 

FIG. 2. The curves display the ratio |[F (k) - F(k)]/F (k) as a function of Z. 
Now let us consider G again with the specific forms of F(k) displayed in Fig. (1-2) 




Note that, given the 1/k 2 dependence of F(k) and of Eq. (39-40), this is a logarithmic 
integral of k (i.e dk/k) times a function of k^/u. Therefore the integration is really over the 
combination variable k^y/ou. Thus falls out of the integral, and the Coulomb correction 
function G does not depend on 7 or u. 
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I have evaluated G numerically and found it exactly converging to Lee and Milstein's 
result according to the expected improved precision with decreasing mesh size. I attained 
agreement to one part in 10 6 . 

Conjecturing that the detail of the cutoff should not matter, I replaced the function 
K (pu/'y) with a different function that also goes as — ha(pcv/2 , y) (plus an irrelevant constant, 
1/2 + ^Euier) for small p and also cuts off exponentially at p ~ 

Upcoh) = ^^[Klif^h) - KoipuhWpu/i)]. (43) 

Calculations of G with L in place of K similarly converge numerically to the result of Lee 
and Milstein with agreement to one part in 10 6 . Note however the non-identical shapes of 
the contribution to G as a function of kj/u for the K Q and L Q transverse potential forms 
exhibited in Fig. (3), even though the area above the two curves (the value of G) is identical. 

Coulomb Corrections 

Region of ky/(0 contributing to G 




-2 



10 



100 



ky/co 



FIG. 3. Region of kj/uj contributing to the Coulomb correction integral G for Z = 82. 
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Now we can begin to understand the result of Lee and Milstein. The reason that "the 
integral G itself is a universal function of Za n is that the first order k 2 factor from the 
expansion makes the integral G logarithmic and so, contrary to what Lee and Milstein 
state, G does not "depend. ..on the regularization parameters (the radius of screening)". 
The radius of screening, i.e. j/co, is finite, but it has fallen out of the problem. Furthermore 
"the main contribution to the integral" does not "come from the region of small k" but, as 
is seen from the plot of the physically motivated K curve in Fig. (3), the main contribution 
is peaked at yk/u = 2.8 and spreads out between half-maxima at 1.3 and 7.5. 

Note that the decoupling of the Coulomb corrections from •y/u seen in G is only valid 
to first order in k. Including higher order terms in k or, alternatively carrying out a full 
numerical evaluation of e.g. Eq. (26), would necessarily restore some dependence on 7/cj 
to the Coulomb corrections. A previous Monte Carlo perturbation theory calculation of 
Bottcher and Strayer [18] displays the pair production cross section as a function of Pt = 

+ q±, and shows a significant deviation between an exact Monte Carlo evaluation of the 
cross section and evaluation using a two peak approximation (in particular see Fig. (9) of 
Ref. [18]). Since in carrying out their calculation, Lee and Milstein made a variety of a two 
peak approximation (assuming Pt = p±_ + qj_ small), one has to assume that the precision 
of their results is limited. 



IV. GENERAL OBSERVATIONS 



To lowest order in transverse momentum (small k and small Pt = P± + q±), Coulomb 
corrections do exist as a universal function of f(aZ), where f(aZ) is the same function of 
Bethe and Maximon derived for Coulomb corrections to electron-positron pair photoproduc- 
tion. These Coulomb corrections reduce the uncorrelated electron or positron production 
cross sections and the number weighted total pair cross section. 

In general and not limited to lowest order in transverse momentum, Coulomb corrections 
are a function of only Z and the combination variable ky/uo. Coulomb corrections arise from 
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the finite cutoff of the transverse spatial integral at and vanish for large kj/iv. 

Since the CERN data cover a large part of the momentum range of produced positrons 
and scale perturbatively, they still seem to present a puzzle. It would be useful to carry 
out full calculations of the total number weighted cross section a T as well as of the uncorre- 
cted momentum dependent electron and positron cross sections da{p) and da(q), utilizing 
the transverse integrals with a correct physical cutoff. Since the CERN data only detects 
positrons, comparison with a full calculation of da(q) is appropriate. 
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